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Abstract. We stress the relevance of the two features of translational invariance 
and atomic nature of the gas in the quantum description of the motion of a 
i-^ ' massive test particle in a gas, corresponding to the original picture of Einstein 

used in the characterization of Brownian motion. The master equation describing 
■J-J ' the reduced dynamics of the test particle is of Lindblad form and complies with 

C ' the requirement of covariance under translations. 

■ 

^ . 

jZj' 1 Introduction 

In recent times there has been growing interest in the study of open quantum systems [l] , driven 
\ by researches on both applications and foundations of quantum theory. Indeed the community 

. of researchers involved in the subject, in spite of the precise label attached to it, ranges from 

physicists and chemists, to mathematicians and probabilists. As it often happens the richness of 
the subject allows for different approaches and biases in the treatment of the same problem. In 
the present paper we come back to the quest for a deep understanding of one of the apparently 
simplest, but also paradigmatic situations in open quantum system theory, that is to say the 
motion of a quantum massive test particle in a homogeneous gas. Two guiding ideas, also present 
in the classical approach to Brownian motion of a suspended particle by Einstein, appear in the 
treatment: symmetry under translations and discrete nature of matter. Following work done 
in [2] we will see that these two concepts reflect themselves in the mathematical properties of 
the mapping describing the reduced dynamics, which has to be covariant under translations, 
as well as in the physical meaning of the two-point correlation function of the gas appearing in 



2 Translational Invariance 

2.1 Microscopic Hamiltonian 

As a first step we characterize the general structure of microscopic Hamiltonians accounting for 
a translationally invariant reduced dynamics for the test particle. We consider a test particle 
subject to a translationally invariant interaction with a homogeneous bath, with a potential 
at most linearly depending on position, e.g., a constant gravitational field. The microscopic 
Hamiltonian may be written as usual 

H = Hs + Hb + V, (1) 
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where the subscripts S and B stand for system and bath respectively, while Hs and Hb satisfy 
the aforementioned constraints. The interaction relies on a translationally invariant potential 
and is of the general form 

V = JdhJ d%As{x)V{x~v)AB[y), (2) 

where As{x) and AB{y) are self-adjoint operators of system and bath respectively, expressing 
the coupling between the two. The invariance under translations of the potential allows to 
express Eq. ([2]) in a very simple way in terms of the Fourier transformed quantities according 
to 

V = j dXlV{Q)As{Q)A^BiQ)- (3) 

In order to focus on a quantum description of Brownian motion we consider a density-density 
coupling, so that As/b{x) = Ng/B{x), with Ns/b{x) the number-density operator of system 
and bath respectively, whose Fourier transform pq 

PQ = j dPxe'T.Q--NB{x) (4) 

is also called Q-component of the number-density operator [3|4j . Eq. ([3]) thus becomes 

V = j d^V{Q)As{Q)plj{Q). (5) 

Note that an interaction of the form ([5]), besides being translationally invariant, commutes 
with the number operators Ns and Nb, so that the elementary interaction events do bring in 
exchanges of momentum between the test particle and the environment, but the number of 
particles or quanta in both systems are independently conserved, thus typically describing an 
interaction in terms of collisions. 



2.2 Quantum linear Boltzmann equation 



The case of density-density coupling given by ([5|), when the reservoir is a free quantum gas, 
has been dealt with in [5l6l7j . and the relevant test particle correlation function turns out to 
be the so-called dynamic structure factor |3|4j 



(6) 



where contrary to the usual conventions, momentum and energy are considered to be positive 
when transferred to the test particle, on which we are now focusing our attention, rather than 
on the macroscopic system. The master equation then takes the form 



dg 



(7) 



=^«-^v/^(Q, E{Q, P))gVS{Q, E{Q, P))e 



{S{Q,EiQ,P)),g} 



where Hq is the free particle Hamiltonian, X and P position and momentum operator for the 
test particle, rigas the density of the homogeneous gas, and the dynamic structure factor appears 
operator-valued: in fact the energy transfer in each collision, which is given by 



(8) 
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with M the mass of the test particle, is turned into an operator by replacing P with P. For 
the case of a free gas of particles obeying Maxwell-Boltzmann statistics the dynamic structure 
factor takes the explicit form 

with /? the inverse temperature and m the mass of the gas particles. 

This result has been confirmed in recent, more general work |8I9| , not relying on the Born 
approximation, where instead of the Fourier transform of the interaction potential the full 
scattering amplitude describing the collisions between test particle and gas particles appear 
operator-valued. In the general case the master equation describing the coUisional dynamics 
takes the following form 

*[Ho,e] (10) 



dt h 



+ Jd^ J [e'^^-^'^'^L (p, P; Q) gL^ (p, P; Q) Q-'Q-^/f^ 

-\{L^ {p,P;Q)Lip,P;Q),g} 
where the Lindblad operators L {p, P; Q) are given by 

^ (P' - \l^f f^^i {p^q^p^q) -^M {p^q.p^q) + (11) 




involving the elastic scattering amplitude / {Pf,Pi) and the momentum distribution function 
(p) of the gas momenta given by the Maxwell-Boltzmann expression. We have further 
denoted relative momenta as rel (p,P) = (m*/m)p— (to*/M)P, with the reduced mass, 
while the subscripts ||Q and ± Q indicate the component of a vector (or operator) parallel and 
perpendicular to the momentum transfer Q, so that P\\q ^ {P ■ Q) Q/Q^ and P±q — P P\\q 
respectively. Exploiting the identity 

m / ^mQ m \ m ( , ( 2mE {Q,P) + Q\ 



= Mmb (p±q) "5" {Q.P\\q) , 

where in the last line /i^^jg {p±q) denotes the Maxwell-Boltzmann distribution over transverse 
momenta, the Lindblad operators can also be written 



L {p, P; Q) = J^f ( rel (p^Q, P^g) - |, rel (p^g, P^q) + | 



Mmb (p±q) Vs (Q,p), 

thus putting again into evidence the appearance of the dynamic structure factor, whose posi- 
tivity has been exploited in order to take the square root. 

The relevant correlation function for the dynamics is thus given by the Fourier transform 
with respect to energy of the time-dependent auto-correlation function of the operator of the 
bath appearing in Eq. ([5]). The appearance of the dynamic structure factor has an important 
physical meaning, linking the dynamics of the test particle to the density fluctuations in the 
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medium, as we shall see in Sect. |3l expressing the molecular, discrete nature of matter, that is 
to say one of the basic insights gained by Einstein's description of Brownian motion. 

The master equation (|10|1 . or ([7|) when considering the Born approximation, can be seen as a 
quantum counterpart of the classical linear Boltzmann equation, as discussed in [2|9] , in that it 
addresses in a quantum framework the same physical situation described by the classical linear 
Boltzmann equation. This is also confirmed by the fact that the diagonal matrix elements in the 
momentum representation of the quantum linear Boltzmann equation do give back the classical 
linear Boltzmann equation, obviously expressed with the quantum scattering cross section. 



2.3 Translation-covariant quantum dynamical semigroups 



In Sect. l2.2l we have considered a test particle interacting through collisions with a homogeneous 
background gas. As stressed in Sect. 12.11 the collisions are to be described by an interaction 
potential only depending on the relative coordinate, so as not to spoil invariance under trans- 
lations. These two requirements lead to a natural general constraint on the structure of the 
mapping giving the reduced dynamics. In fact homogeneity of the bath implies that the statis- 
tical operator describing its equilibrium state commutes with the momentum operator of the 
bath. Similarly the considered interaction V ensures that the total Hamiltonian H commutes 
with the momentum operator of the whole system, which we can write as Ps + Pb- Let us 
now consider the reduced operator of the test particle at time t obtained by taking the trace 
over the bath degrees of freedom of the statistical operator of the total system. Considering a 
factorized initial state one has 

Ut [gs] - Trs [e-^"'gs ® ^56+^^*) . 

Exploiting further the aforementioned constraints 



[gB,PB]=0 and [H,Ps + Pb]=0, 
one immediately has, for any vector a £ R"^ 



so that a mapping Ut giving the reduced dynamics must obey 



This condition is known as covariance under translations. Focusing on the Hilbert space of the 
massive test particle considered in the present paper, that is to say {^'^) 7 the condition 
can be stated as follows. Given the unitary representation U(a) — exp(— ia • P/h), a £ R'^ 
of the group of translations IR'^ in '"^ mapping C acting on the statistical operators 

in this space is said to be translation-covariant if it commutes with the action of the unitary 
representation, i.e. 

/:[U(a)gUt(a)] = U(a)/:[g]U^(a), (13) 

for any statistical operator g and any translation a. The general structure of generators of 
quantum dynamical semigroups complying with this covariance condition has been obtained 
by Holevo |10lll| . and it turns out that the requirement of translation covariance puts very 
stringent constraints on the Lindblad operators appearing in the expression of the generator. 
These results, while obviously fitting in the general framework set by the famous Lindblad 
result [12|13j . go beyond it giving much more detailed information on the possible choice of 
operators appearing in the Lindblad form, information conveyed by the symmetry requirements 
and relying on a quantum generalization of the Levy-Khintchine formula. They therefore also 
provide a precious starting point for phenomenological approaches exploiting relevant physical 
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symmetries. Referring to the papers by Holevo for the related mathematical details (see also 
[l4] for a brief resume) , the generator can be expressed as 



C[g] = --[H{P),e]+CG[Q\+Cp[Ql 



(14) 



with H{P) a self-adjoint operator which is only a function of the momentum of the test particle. 
The so-called Gaussian part Cg is given by 



Yo + i7eff(X,P),e]+^ 



k=l 



l{KlKk,g} 



where 



(15) 



Kk = Yfe + Lfe(P), Yfc = P) = 2z - 4(P)Yfc) 



1=1 fc=i 
with fc = 0, . . . , r < 3 and Uki G R, while the remaining Poisson part takes the form 

Cp[g] = I df^iQ) J2 e^'^-^i:,(Q, P)gL]{Q, P)e-^'3 >^ - - [l]{Q, P)L,{Q, P), g} 



(16) 



with dfi(Q) a positive measure. The names Gaussian and Poisson arise in connection with the 
different contributions in the classical Levy-Khintchine formula [13]. In the Gaussian part the 
Yfe are linear combinations of the three position operators of the test particle, while the generally 
complex functions Lk{P) have an imaginary part accounting for friction, typically given by a 
linear contribution, corresponding to a friction term proportional to velocity. In the Poisson 
part a continuous index Q appears, together with the usual sum over a discrete index j. The 
expression is characterized by the appearance of the unitary operators exp(iQ-X/fi,), expressing 
momentum kicks, and of the functions Lj{Q, P), operator- valued in that they depend on the 
momentum operators of the test particle P. 



3 Fluctuation-dissipation theorem 

As already stressed the two-point correlation function appearing operator-valued in the master 
equation is the dynamic structure factor ([6]) , where the Fourier transform of the number-density 
operator pQ, as given in (j4|), appears. This function is directly related to the density fluctuations 
in the medium, as it can be seen writing it in the following way ^3J: 

S{Q,E) = ±^j dtj d^r e*(^*-Q-)l j d\ {NB{y)NB{x + y,t)), (17) 

i.e., as Fourier transform with respect to energy and momentum transfer of the time depen- 
dent density correlation function . Here the connection with density fluctuations and therefore 
discrete nature of matter is manifest. Introducing the real correlation functions 

r{Q,t)^J^{[PQ{t),ph]) and 4,+ {Q,t)^^{{pQ{t),pl^}), (18) 

where {, } denotes the anticommutator, the fluctuation-dissipation theorem can be formulated 
in terms of the dynamic structure factor as follows 

2 

</,-(Q,t) = -- / dE sm{Et/h){l~e'^^) S{Q,E) (19) 

J — oo 

2 f° 

0+(Q, t) = -- dE cos {Et/h) coth {I3/2E) (l - e'^^) S{Q, E). 
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We recall that contrary to the usual perspective in linear response theory we take as positive 
momentum and energy transferred to the particle. The dynamic structure factor can also be 
directly related to the dynamic response function x"(Qi E) f4j, according to 



1 - coth ( 



^"iQ,E) = --^x"iQ,E), (20) 



the relationship leading to the important fact that while the dynamic response function is 
an odd function of energy, the dynamic structure factor obeys the so-called detailed balance 
condition 

S{Q,E) = e-'^'^S{-Q,-E), (21) 

a property granting the existence of a stationary state for the master equation [7j . 

The significance of the appearance of the dynamic structure factor in connection to the 
so-called fluctuation-dissipation theorem is to be traced back to a seminal paper by van Hove 
[161 17| . In fact he showed that the scattering cross-section of a microscopic probe off a macro- 
scopic sample can be written in Born approximation 

.iP) = {2nhf(^X^\Vm'S{Q,E), (22) 



dfip,dEpr ' ' ' yiTxh^ ) p 

where a particle of mass M changes its momentum from P to P' = P + Q scattering off 
a medium with dynamic structure factor S{Q,E). This can be seen as a formulation of the 
fluctuation-dissipation relationship for the case of a test particle interacting through collisions 
with a macroscopic fluid. The energy and momentum transfer to the particle, characterized 
by the expression of the scattering cross-section at l.h.s. of (I22p are related to the density 
fluctuations of the macroscopic fluid appearing through the dynamic structure factor at r.h.s. 
of ((22|) . One of the basic ideas of Einstein's Brownian motion, i.e., the discrete nature of matter, 
once again appears in the formulation (|22p of the fluctuation-dissipation relationship. 



4 Friction coefficient for quantum description of Brownian motion 

We now come to the master equation for the quantum description of Einstein's Brownian 
motion. The requirement of translational invariancc has been settled in Sect. [21 while the con- 
nection between reduced dynamics of the test particle and density fluctuations in the medium, 
coming about because of its discrete nature, has been taken into account in Sect. [3] The last 
step to be taken is to consider the test particle much more massive than the particles making up 
the gas, i.e., the Brownian limit m/M <^ 1, which in turn implies considering both small energy 
and momentum transfers, similarly to the classical case |18|7j . We therefore start from ^ and 
consider a free gas of Maxwell-Boltzmann particles, so that taking the limiting expression of © 
when the ratio between the masses is much smaller than one leads, of necessity as can be seen 
from the Gaussian contribution in Holevo's result (fTS]) but also from previous work [19|20|21] . 
to a Caldeira Leggett type master equation, however without shortcomings related to the lack 
of preservation of positivity of the statistical operator. The master equation takes the form 

f - [Ho, [X. {P. [X. - ^ E [P^, [P^, 9\] , (23) 

i—l i—1 i—1 

with 

M ph^ 

Dm = -rrV and Drr = V- (24) 

The friction coefficient r] is uniquely determined on the basis of the microscopic information on 
interaction potential and correlation function of the macroscopic system, according to 

^=^y(2^^)Xas/ dXi\Vm'^S{Q,E = 0), (25) 
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the factor 3 being related to the space dimensions, thus proving in a specific physical case of 
interest the so-called standard wisdom expecting the decoherence and dissipation rate to be 
connected with the value at zero energy of some suitable spectral function 22J. The transition 
from the general master equation ([7]) to the approximate expression ()23|) has been considered 
in detail in [23] . where the microscopic expression for the friction coefficient has been worked 
out in detail for the case of a constant scattering cross section. In order to point out the 
connection with classical Brownian motion as described by Einstein we stress that ()23|) is 
a quantum version of the classical Kramer's equation. This can be easily seen considering 
the usual correspondence rules between classical and quantum mechanics, sending position to 
multiplication by the variable and momentum to derivation, or also considering the expression 
of for the Wigner function, which reads 

^wix,p) = -^-yxW{x,p) 

+7?Vp • {PW (X, P)) + DppApW {X, P) + D^^AxW {X, P) , 

and in the strong friction limit leads to the classical Smoluchowski equation with a small 
quantum correction [24 j . 

Acknowledgements 

The work was partially supported by the Itahan MIUR under PRIN05 (BV). 



References 

1. H.P. Breuer, F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press, 
Oxford, 2007) 

2. F. Petruccione, B. Vacchini, Phys. Rev. E 71, 046134 (2005) 

3. S. Lovesey, Theory of neutron scattering from condensed matter. Vol.1. Nuclear scattering (Claren- 
don Press, Oxford, UK, 1984) 

4. L. Pitaevskii, S. Stringari, Bose-Einstein condensation (Oxford University Press, Oxford, 2003) 

5. B. Vacchini, Phys. Rev. Lett. 84, 1374 (2000) 

6. B. Vacchini, Phys. Rev. E 63, 066115 (2001) 

7. B. Vacchini, J. Math. Phys. 42, 4291 (2001) 

8. K. Hornberger, Phys. Rev. Lett. 97, 060601 (2006) 

9. K. Hornberger, B. Vacchini (2008), to appear in Phys. Rev. A [arXiv:quant-ph/0711.3109] 

10. A.S. Holevo, Rep. Math. Phys. 32(2), 211 (1993) 

11. A.S. Holevo, Rep. Math. Phys. 33(1-2), 95 (1993) 

12. G. Lindblad, Comm. Math. Phys. 48(2), 119 (1976) 

13. V. Gorini, A. Kossakowski, E.C.G. Sudarshan, J. Math. Phys. 17(5), 821 (1976) 

14. B. Vacchini, Int. J. Theor. Phys. 44(7), 1011 (2005) 

15. W. Feller, An introduction to probability theory and its applications. Vol. II (John Wiley & Sons 
Inc., New York, 1971) 

16. L. Van Hove, Phys. Rev. 95, 249 (1954) 

17. F. Schwabl, Advanced quantum mechanics, 2nd edn. (Springer, New York, 2003) 

18. C.S.W. Chang, G.E. Uhlenbeck, The kinetic theory of gases, in Studies in statistical mechanics, 
edited by J.D. Boer (North-Holland, Amsterdam, 1970), Vol. 5 

19. G. Lindblad, Rep. Mat. Phys. 10(3), 393 (1976) 

20. A. Sndulescu, H. Scutaru, Ann. Physics 173(2), 277 (1987) 

21. L. Diosi, Europhys. Lett. 30, 63 (1995) 

22. R. Alicki, Open Syst. Inf. Dyn. 11(1), 53 (2004) 

23. B. Vacchini, K. Hornberger, Eur. Phys. J. ST 151, 59 (2007) 

24. B. Vacchini, Phys. Rev. E 66(2), 027107 (2002) 



